Atmospheric scintillation at Dome C, Antarctica: implications for 

photometry and astrometry 



O 
O 
(N 

;h 

< 
(N 



> 

oo 
m 
in 

o 

o 

o 



1 



S.L. Kenyon, J.S. Lawrence, M.C.B. Ashley, J.W.V. Storey 

School of Physics, University of New South Wales, Sydney 2052, Australia 

suzanneSphys . unsw . edu . au 

A. Tokovinin 
Cerro-Tololo Inter American Observatory, Casilla 603, La Serena, Chile 

and 

E. Fossat 

Lahoratoire Universitaire d'Astrophysique de Nice, Universite de Nice, France 
Submitted to PASP 6th April 2006, accepted 26th April 2006 



ABSTRACT 

We present low- resolution turbulence profiles of the atmosphere above Dome C, Antarctica, 
measured with the MASS instrument during 25 nights in March - May 2004. Except for the 
lowest layer, Dome C has significantly less turbulence than Cerro Tololo and Cerro Pachon. In 
particular, the integrated turbulence at 16 km is always less than the median values at the two 
Chilean sites. From these profiles we evaluate the photometric noise produced by scintillation, 
and the atmospheric contribution to the error budget in narrow-angle differential astrometry. In 
comparison with the two mid-latitude sites in Chile, Dome C offers a potential gain of about 
3.6 in both photometric precision (for long integrations) and narrow-angle astrometry precision. 
These gain estimates are preliminary, being computed with average wind-speed profiles, but the 
validity of our approach is confirmed by independent data. Although the data from Dome C 
cover a fairly limited time frame, they lend strong support to expectations that Dome C will 
offer significant advantages for photometric and astrometric studies. 

Subject headings: Astrometry — atmospheric effects — site testing — turbulence 



1. Introduction 

The potential of the Antarctic plateau for as- 
tronomy has been recognized for many years. In 
the early 1990's Gillingham (1991, 1993) sug- 
gested that the atmospheric turbulence above a 
thin boundary layer would generally be very weak. 
This was confirmed by the first measurements at 
the South Pole by Marks et al. (1999). These 
showed that most turbulence is found close to the 
surface; in this case confined to a 200 m thick 



layer, with very weak turbulence at high altitudes. 
In comparison, at most temperate sites the tur- 
bulence is strong in the tropopause and above, 
caused by the interactions of the jet stream with 
temperature gradients in the tropopause. A brief 
history of astrophysics in Antarctica is presented 
by Indermuehle et al. (2005). 

Dome C, Antarctica is potentially one of the 
best astronomical sites in the world. As a local 
maximum in elevation on the plateau. Dome C 
enjoys very low surface wind speeds, on average 



2.9 m s-i(Aristidi et al. 2005). 

The atmospheric turbulence at Dome C has 
now been studied with four different techniques: 
acoustic radar, MASS (a scintillation profiling 
technique; see later), DIMM, and microthermal 
sensors. An acoustic radar, or SODAR, emits 
sound pulses into the air and derives the strength 
of the atmospheric turbulence as a function of 
height from the intensity and delay time of the 
reflected sound. A DIMM, or Differential Image- 
Motion Monitor, observes the relative motion of 
two images of the same star viewed through two 
sub-apertures of a small telescope. From this, 
the DIMM can derive the integrated atmospheric 
seeing. Microthermal sensors, carried aloft on a 
weather balloon, make an in situ measurement of 
the temperature fluctuations of the air as a func- 
tion of height all the way to the top of the atmo- 
sphere. 

SODAR measurements in the early months of 
2003 by Travouillon (2005) showed that, as ex- 
pected, the surface turbulent layer at Dome C 
was much thinner that at South Pole. Combined 
MASS and SODAR measurements of the turbu- 
lence in winter 2004 gave an average seeing of 
0.27" above 30 m, with the seeing below 0.15" for 
25% of the time (Lawrence et al. 2004a). DIMM 
measurements in winter 2005 (Agabi et al. 2006) 
confirmed these results, showing an average see- 
ing of 0.25" above the ground layer (Agabi et al. 
2006, Fig. le), and balloon microthermal mea- 
surements by the same authors implies a median 
seeing of 0.36" ± 0.19" at a height of 30 m. They 
also showed the existence of an intense turbulent 
boundary layer, finding a median seeing of 1.9" 
from ground level, with 87% of the total atmo- 
spheric turbulence confined to the first 36 m of 
atmosphere. 

As expected, at Dome C the sky background 
in the infrared (Walden et al. 2005) is lower than 
at temperate sites because of the extremely cold 
temperatures and lower precipitable water vapor. 
For the same reason, the atmospheric transmission 
in the sub-millimeter is higher than at temperate 
sites (Calisse et al. 2004). An assessment of the 
optical sky brightness has been recently published 
by Kenyon & Storey (2006). 

To date, the winter time scintillation at Dome C 
has only been estimated from atmospheric models 
(e.g., Swain et al. 2003). Scintillation is an impor- 



tant factor in measurements requiring high pre- 
cision photometry (e.g., extra-solar planet detec- 
tion) and astrometry, and of objects with very fast 
intensity changes (e.g. astroseismology; Heasley 
et al. 1996; Fossat 2005). 

Here we evaluate the scintillation noise contri- 
bution to photometry and the atmospheric noise 
contribution to narrow-angle astrometry, using a 
set of low-resolution turbulence profiles measured 
at Dome C in 2004. The instrument and data are 
described in § 2. In § 3 we outline the theory of 
atmospheric turbulence, scintillation and interfer- 
ometry. In § 4 we present our results and discus- 
sion. 

2. MASS Measurements 

The turbulence profile of the atmosphere at 
Dome C was monitored with a MASS (Multi- 
Aperture Scintillation Sensor) instrument during 
the first 2 months of the 2004 night time: 23 March 
2004 to 16 May 2004. The analysis of these data 
in terms of seeing has been reported by Lawrence 
et al. (2004a). 

The MASS instrument and theory are described 
in detail in Kornilov et al. (2003) and Tokovinin 
et al. (2003). In brief, starlight is directed via a 
telescope onto four concentric annular mirrors that 
split the entrance aperture into rings with pro- 
jected outer diameters of 19, 32, 56, and 80 mm. 
Each of the four beams is directed to a minia- 
ture Hamamatsu photomultiplier which samples 
the stellar intensity at a 1 kHz rate. Four normal 
and six differential scintillation indices are calcu- 
lated for each 1 second integration and further av- 
eraged during 1 minute. The set of 10 indices is 
fitted to a model of six fixed layers at heights 0.5, 
1, 2, 4, 8 and 16 km above the observatory. For 
each layer i, the integrated turbulence Ji (m^/'^) 
is calculated: 



J,.= 



Layer i 



Cl{h)dh, 



(1) 



where C^{h) (m~^/^) is the refractive index struc- 
ture constant and h is the height above the site. 
The spectral response of MASS is from 400 to 
550 nm with a FWHM bandwidth of 100 nm. 

The profile restoration from scintillation indices 
is a delicate procedure, and errors may reach 10% 
of the total turbulence integral. The errors are 



larger for the lower layers, while the 2 highest lay- 
ers (8 km and 16 km) are measured well. The 
second moment of the turbulence profile used in 
this paper is measured by MASS with high reli- 
ability, as well as lower moments. This has been 
demonstrated by inter-comparing MASS and SCI- 
DAR instruments (Tokovinin ct al. 2005). 

The Dome C MASS (Lawrence et al. 2004b), 
operated in the AASTINO (Automated Astro- 
physical Site Testing International Observatory; 
Lawrence et al. 2005), uses a gimbal-mounted 
siderostat mirror feeding a fixed 85 mm refract- 
ing telescope. 

MASS instruments also operate at a number 
of other sites. To provide a comparison to the 
Dome C results, we have also included the publicly 
available data^ from the Cerro Tololo and Ccrro 
Pachon observatories in Chile. The profiles for 
Cerro Pachon have been discussed and modeled 
by Tokovinin & Travouillon (2006). 

3. Theory 

3.1. Turbulence 

Many astronomical measurements are limited 
by the Earth's atmosphere. A wavefront located 
at height h and horizontal position vector x in 
the atmosphere can be described by its complex 
amplitude ^^(x) (Roddier 1981), 



*h(x) = c^ 



= pXh(x)+ii/'h(x) 



(2) 



where Xh (x) is the logarithm of the amplitude and 
V'/i(x) is the phase of the wave. 

Atmospheric turbulence introduces pure phase 
distortions. As the wavefront propagates through 
the atmosphere, amplitude modulations appear as 
well. In the geometric optics approximation, phase 
perturbations act as positive or negative lenses, 
changing the wavefront curvature and producing 
intensity modulation at the ground. Diffraction 
is also important and defines the size of the most 
effective atmospheric "lenses" to be of the order 
of the Fresnel radius, 



rp ~ (Xh) 



1/2 



(3) 



1 Obtained from the National Optical Astronomy 
Observatory web page "Sites Data Access," at 
http://139.229.ll.21/ 



where A is the wavelength of light and h is the 
height of the turbulent layer above the observatory 
site. For example, if the dominant turbulence layer 
is at 10 km then at 500 nm, rp ^ 7 cm. 

Dravins et al. (1997a, b, 1998) present detailed 
discussions of stellar scintillation, including sta- 
tistical distributions and temporal properties, de- 
pendence on wavelength and effects for different 
telescope apertures. 

3.2. Scintillation noise 

The scintillation index aj is used as a mea- 
sure of the amount of scintillation and is defined 
(for small intensity fluctuations) as the variance of 
A J/ (/). In the weak-scintillation regime, af ^ 1, 
the effects of all turbulence layers are additive. 
In this case, the scintillation index is related to 
the refractive index structure constant C^{h) by 
(Krause-Polstorff et al. 1993; Roddier 1981) 



Cl{h)W{h)dh, 



(4) 



where the weighting function W{h) is given by 



W{h) ^IGtt^O.OSS 



x/ \Aif)\\r^/^sin'f^]df. 



(5) 



Here, h is the height above the observatory, A 
is the wavelength, / is the spatial frequency and 
|^(/)| is an aperture filtering function. This ex- 
pression is valid for monochromatic light and has 
to be modified for wide-band radiation. 

For telescope apertures with diameter D <^rF, 
the monochromatic scintillation index is (Roddier 
1981) 

a] = 19.2A-^/6(cos7)""/' / h'>"'Cl{h)Ah, 

Jo 

(6) 

where 7 is the zenith angle. 

The scintillation index for a large circular aper- 
ture with diameter D ^ rp is: 

/•oo 

a| = 172?-^/3(cos7)-3 / h^Ct{h)dh. (7) 



Large apertures effectively average small-scale 
intensity fluctuations, so that only atmospheric 



lenses of the order of the aperture diameter D 
contribute to the flux modulation. In this case, 
geometric optics applies and the scintillation be- 
comes independent of both the wavelength and 
the spectral bandwidth. 

The above expressions are for very short 
time scale exposures. For exposure times that 
are longer than the time taken for a scintilla- 
tion pattern to cross the telescope aperture (i.e 
t > (TrD)/V±, where D the telescope diameter and 
V± the speed of the turbulence layer) , the scintil- 
lation index can be calculated from (Dravins et al. 
1998) 



aj{t) 



P{h')smc (Tnytjdh', 



(8) 



where v (s^^) is the temporal frequency and P{iy) 
is the temporal power spectrum, given by Yura & 
McKinlcy (1983, and ref. therein) as 

P(^) « 8.27fc2/3 r ^M^Q(^h)dh (9) 
Jo V±{h) 

at the zenith, where k = 2tt/X and 



Qih) 



\Ai.U,fy)\' 



<{h) 



-11/6 



sm 





(10) 




da;. 



Here, /, = i^lV^{h), fy = {2k/h)^/^x, ^h) = 
{2k/hY'^V^{h) and Aif) = {2MfD/2])/{fD/2) 
for a circular aperture with diameter D and /^ = 

J X I J y 

For large t, Dravins et al. (1998) simplify Equa- 
tion 8 to 



aUt) 



m 

2t ' 



(11) 



In the limit of large apertures D ^ rp, we can 
replace the sine in Equation 10 by its argument. 
By setting ly ^ and introducing a new variable 
y ~ fyD, we can show that 

P(0) « 2l.iD-^'^ r ^"^^j^' dfe (12) 

Jo Kl(/i) 

and is independent of the wavelength. 



For a particular set of turbulence and wind pro- 
files, the scintillation noise ai at zenith can be ex- 
pressed as 



(Ji 



where 




D <^rF 

D:§>rF,t:$> inD)/Vi_ 
(13) 



Si - 



S. = 



19.2A-^/6 / h^/^cl{h)dh 



1/2 



17.3 



10.7 



h^Cl{h)dh 



V±(h) 



1/2 



1/2 



dh 



,(14) 
(15) 
(16) 



The scintillation error can be expressed in magni- 
tudes as a/(mag) = 2.51og(cr7 -I- 1). 

In all cases the scintillation noise is dominated 
by the high-altitude turbulence, more so in the 
case of large apertures because of the h^ weighting. 
It is the large-aperture case that is generally of 
more relevance to astronomical photometry. 

3.3. Astrometric interferometry 

The Antarctic plateau has been recognized as 
a potentially favorable site for interferometry be- 
cause the high-altitude turbulence is very weak 
(Lloyd et al. 2002). In particular, high precision, 
very-narrow angle differential astrometry should 
be attainable at Dome C using long baseline in- 
terferometry techniques. This would benefit a 
number of science programs, including extra-solar 
planet searches and the study of close binary and 
multiple star systems (for other examples, see 
Swain et al. 2003; Lloyd et al. 2002 and Sozzetti 
2005). 

Differential astrometric measurement requires 
simultaneous observations of the target and refer- 
ence object. To achieve this, each telescope has a 
dual feed to direct the beam from each star to the 
beam combiner (Shao & Colavita 1992). On com- 
bination of the beams, a fringe pattern is produced 
if the difference between the optical path lengths 
from each arm of the interferometer to the beam 
combiner is within A^/AA (Lane & Muterspaugh 
2004) . The difference between the fringe positions 
of the two stars is measured. Phase referencing 



can be used to improve the limiting magnitude of 
the interferometer if the target star and reference 
object are within the isoplanatic patch (Shao & 
Colavita 1992). 

Uncertainties in astrometric position measure- 
ments arise from instrumental effects (noise, sys- 
tematic) and atmospheric effects associated with 
temporal incoherence and anisoplanatism. See 
Shao & Colavita (1992); Sozzetti (2005) and Lane 
& Muterspaugh (2004) for further details. 

The variance in an astrometric position mea- 
surement caused by anisoplanatism (assuming a 
Kolmogorov turbulence spectrum) is described by 
Shao & Colavita (1992) as 



5.25t 



B-^/^l 




oo Cl{h)h^ 
V{h) 







V(h) 



dh 



Ah Case 1 
Case 2, 

(17) 
where t is the integration time, is the angular 
separation between two stars, C^{h) and V{h) are 
the vertical turbulence and wind profiles, h is the 
height above the site and B is the baseline or di- 
ameter of the entrance pupil. These formulae are 
only approximate, but the exact coefficient is not 
needed for the purpose of site intcr-comparison. 

Case 1 applies to interfcrometry when the in- 
tegration time t ^ B/V and 9h <C B; where h 
and V are the turbulence-weighted effective atmo- 
spheric height and wind speed. Because of the h^ 
weighting, (Tatm in this regime is very sensitive to 
high-altitude turbulence. 

Case 2 is applicable to single dish astrometry 
and is independent of the size of the telescope 
when 9h:§> B and t > 9h/V. 

For a particular set of turbulence and wind pro- 
files, the error CTatm (arcseconds) can be expressed 
as 



CatE 



Cit-1/205-2/3 Case 1 

C2t"i/2g|i/3 Case 2, 



where 



Ci = 472 000 



V{h) 



(18) 



(19) 



and 



C2 = 472 000 



uo 



V{h) 



Ah 



1/2 



.(20) 



Note that the expression for Ci contains the same 
combination of atmospheric parameters as the ex- 
pression for the photometric error 6*3. This is not 
a coincidence, as both narrow-angle astrometry 
and large-aperture photometry are affected by the 
same physical phenomenon - large-scale curvature 
fluctuations of wavefronts. Hence, scintillation in 
large apertures contains information on the po- 
tential accuracy of narrow-angle astrometry at a 
given site. 

4. Results 

Using the eight weeks of MASS data from 
Dome C, we extracted 11532 turbulence profiles 
spread over 51 nights (we use "night" to mean the 
period within 24 hours when the Sun is further 
than 10° below the horizon). These data were fil- 
tered according to the criteria: Bd/Fjj < 0.03, 
Fd > Fiimit, SFd < 0.003, 0.7 < Fc/Fd < 0.9 
and x^ < 100, where Fjj and Bd are the star 
flux and background measurements in aperture D 
(largest aperture). -Fiimit is a cut-off flux limit set 
to 100 counts for Alpha Trianguli, and 200 counts 
for Beta Crucis and Beta Carinae. The parame- 
ter SF]j shows slow flux variations, used here to 
eliminate data affected by the guiding errors. The 
flux ratio Fc/Fd serves to control the aperture vi- 
gnetting by the entrance window, which was some- 
times covered by snow or frost. The x^ is a mea- 
sure of the fit quality. After filtering, 1853 profiles 
over 26 nights remained for further analysis. 

Each profile includes the integrated turbulence 
Ji (see Equation 1) in layers centered at eleva- 
tions of 0.5, 1, 2, 4, 8 and 16 km above the site 
with vertical resolution Ah/h ^ 0.5. Finally, we 
calculated the scintillation noise and implied as- 
trometric error from each profile. In this section, 
we compare these results with similar data for the 
Cerro Tololo and Cerro Pachon observatories in 
Chile (see Table 1 for information on each data 
set). 

4.1. Turbulence profiles 

Figure 1 shows the cumulative probability that 
the integrated turbulence for each height is less 
than the given Ji. 

Cerro Tololo has the lowest turbulence in the 
0.5 km layer. At 1 km the turbulence at Dome C is 
so low that, for most of the time, it cannot be reli- 



Table 1 
Data sets 



Site 


Location 


Altitude 


Date 


range 


Number of 
Nights Profiles 


Dome C 
Cerro Tololo 
Cerro Pachon 


123°21' E 75°6' S 
70°48' W 30°9' S 
70°44' W 30° 14' S 


3260 m 
2215 m 
2738 m 


23 Mar 2005 - 
19 Mar 2002 
9 Jan 2003 - 


- 16 May 2004 

- 2 Feb 2006 

- 30 Jan 2006 


26 
573 
293 


1853 
98887 
39819 



Note. — We use "night" to mean the period within 24 hours when the Sun is further than 10° below 
the horizon 



ably measured with MASS. Dome C has a slightly 
higher probability of smaller turbulence in the 4 
and 8 km layers. However, the most significant 
difference between the sites is in the 16 km layer; 
at Dome C the integrated turbulence in this high- 
altitude layer is always less than the median values 
at Cerro Tololo and Cerro Pachon. 

The Cerro Tololo and Cerro Pachon sites are 
only 10 km apart and have a 400 m altitude dif- 
ference. Hence, we expect identical high-altitude 
turbulence for those sites. The differences seen 
in Figure 1 refiect mostly different seasonal cover- 
age of the data sets (more winter-spring data for 
Cerro Tololo) coupled to the systematic seasonal 
trends in high-altitude turbulence. Similar cau- 
tion is warranted for the Dome C data that cover 
only 25 nights. 

4.2. Wind profiles 

The photometric error for long integration 
times (equations 13 and 16) and the astromet- 
ric errors (equations 19 and 20) depend not only 
on the turbulence but also on the wind speed 
profile. As these quantities are likely correlated, 
the correct way to estimate the errors requires 
simultaneous data on wind and turbulence. The 
wind profiles can, in principle, be retrieved from 
the global meteorological databases like NCEP 
(National Centers for Environmental Prediction). 
However, here we adopt a simplified approach and 
use fixed wind profile models instead. Hence, the 
distributions derived here may be not realistic. 

Owing to the strong h^ weighting, the astro- 
metric and photometric errors arc almost entirely 



determined by the highest MASS layer at 16 km. 
Hence, the adopted wind speed in this layer criti- 
cally influences our results. 

The wind speed profiles for Cerro Pachon 
and Cerro Tololo were modeled using a constant 
ground layer speed Vg plus a Gaussian function to 
represent the jet stream contribution (Greenwood 
1977). 



V{h) = Vg + Vt exp 



h-H 
T 



(21) 



where h is the altitude above the observatory. We 



set V„ = 8 



Ft = 30 m s-\ if = 8 km and 



T ~ 4: km by comparing the model to the Cerro 
Pachon wind profiles in Avila et al. (2000) and 
Avilaet al. (2001). 

The summer wind speed profile at Dome C also 
shows a Gaussian peak at the (somewhat lower) 
tropopause layer {'^ 5 km) and fairly constant 
wind speed at other elevations (see Figure 4 of 
Aristidi et al. 2005). In the winter, the wind 
speed profile is different, showing an increase in 
stratospheric wind speeds and no peak at the 
tropopause. So far only 3 profiles of the win- 
ter wind speed have been published (Agabi et al. 
2006) . Figure 2 shows the average winter and sum- 
mer time wind profiles, we used the winter-time 
wind profile in this work. 

4.3. Scintillation noise 

The scintillation noise was calculated for the 
three regimes discussed in § 3.2, using the results 
from each site. Figure 3 shows the cumulative 
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Fig. 1. — Cumulative probabilities that the inte- 
grated turbulence for each height (above the sur- 
face) is less than the given Ji for Dome C (solid), 
Cerro Tololo (dashed) and Cerro Pachon (dotted). 
The large fraction of very low Ji values for the 
0.5 km and 1 km layers is an artifact of the MASS 
profile restoration method in situations when these 
layers do not dominate. 




10 20 30 

Wind Speed (m/s) 

Fig. 2. — Average wind speeds profiles at Dome C 
in the winter (Agabi et al. 2006) and summer 
(Aristidi et al. 2005), and a model of the wind 
profile at Cerro Tololo and Cerro Pachon. 

probabilities for 5*1, 52 and 5*3. 

For short time scales, the median scintillation 
noise at Dome C is a factor of ~ 2 less than 
at Cerro Tololo and Cerro Pachon, in the small 
aperture regime. For larger apertures the gain is 
slightly higher, ~ 2.4, because of the weaker high 
altitude turbulence at Dome C. As an example, 
for a 4 m diameter telescope, the median values of 
the scintillation noise at each site are: 1.2 mmag 
(Dome C), 3.2 mmag (Cerro Tololo) and 2.8 mmag 
(Cerro Pachon). 

The more relevant figure is the scintillation 
noise for long exposure times. We used the wind 
speed models discussed above, to calculate this pa- 
rameter. Dome C offers a potential gain of about 
3.6 in photometric precision compared to Cerro 
Tololo and Cerro Pachon. From the results we 
calculate the median photometric error expected 
on a 4 m telescope for i = 60 s to be: ~ 53 /imag, 
at Dome C; ^ 200 /imag at Cerro Tololo, and 
^ 180 /imag at Cerro Pachon. 

As a comparison, Dravins et al. (1998) mea- 



sured P(j^) at La Palma using various small aper- 
tures. From their results they extrapolate P{0) — 
5 X 10~^ s for a 4 m aperture at zenith, which gives 
(7/ = 220 ^mag for a 60 s integration, similar to the 
typical values for Cerro Tololo and Cerro Pachon. 
Our results are also consistent with those mea- 
sured at Kitt Peak and Mauna Kea by Gilliland 
et al. (1993). 

4.4. Astrometry 

The constants Ci and C2 (§3.3) were calculated 
for each site; cumulative probabilities are shown 
in Figure 4. The median astrometric error (Tatm at 
Dome C is ~ 3.5 times less than the median values 
at Cerro Tololo and Cerro Pachon. In Figure 5, 
c^atm at Dome C is plotted for several baselines, as 
a function of separation angle 9 for an integration 
time of 1 h. 

We note that the advantage of Dome C for 
narrow-angle astrometry over mid-latitude sites is 
even larger than its advantage in the fast scintilla- 
tion. This difference is related to the adopted wind 
speed at 16 km altitude (27 m s^^ and 8.5 m s^^ 
for Dome C and Cerro Pachon, respectively). Tur- 
bulence at Dome C is known to be slow (large time 
constant, see below), but the h^ weighting in the 
expressions for the photometric and astrometric 
errors reverses this conclusion because the high 
altitude turbulence dominates the calculation. 

Our conclusions are conditional on the adopted 
wind-speed models. Using the mean Dome C win- 
ter wind speed we calculated a median C'l value of 
140 arcsec rad~^m^/'^s^/^; decreasing the 16 km 
wind speed to 7 m s~^ gave a median Ci value 
of 200 arcsec rad^^m^/'^s^/^, still well below the 
median values at Cerro Tololo and Cerro Pachon. 
As an additional check, we computed Ci from a 
set of six balloon profiles of C^ and wind mea- 
sured at Cerro Pachon in October 1998 (see Avila 
et al. 2000 and Avila et al. 2001 for the discus- 
sion of these data). The Ci values range from 
380 to 660 arcsec rad~^m^s^/^, with a median 
of 480. This is close to the median value for 
Cerro Pachon given in Figure 4. Shao & Colavita 
(1992) calculate Ci at Mauna Kea to be 300 arc- 
sec rad~^m^s^/^, using the results from 2 short 
observing campaigns. 

The fringe phase of an interfcrometric measure- 
ment must be determined within the atmospheric 
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Fig. 3. — Cumulative probabilities of the con- 
stants (a) S'l, (b) 52 (mV6) and (c) ^3 (m^/s s^/^) 
for Dome C, Cerro Tololo and Cerro Pachon. The 
scintillation noise aj is the standard deviation of 
A/// where / is the stellar flux, and is equal 
to: Si for D < rp; 820-^^ for D > rp, and 
SaD-'^/H-^/^ for £> > rj. and t > i7rD)/VjL, 
where D is the telescope diameter, rp is radius of 
the Fresnal zone and t is the integration time. 



coherence time. Table 2 shows the median coher- 
ence times To (Roddier et al. 1982a) and isopla- 
natic angles 9o (Roddier et al. 1982b) at each site 
for wavelengths 500 nni and 2.2 fiia. The me- 
dian coherence time at Dome C, measured with 
MASS, is 7.2 ms at A = 500 nm and 42 ms at 
A = 2.2 /im. Using the measured turbulence pro- 
files and assumed wind profile at Dome C we cal- 
culated To — 9.4 ms at A = 500 nm. 

Phase referencing during the measurement 
(Shao & Colavita 1992) effectively increases the 
coherence time, with the condition that the target 
and reference objects are within the same iso- 
planatic patch. The median isoplanatic angle at 
Dome C is '--^ 3 times larger than at Cerro Tololo 
and Cerro Pachon, allowing wider fields to be used 
for phase referencing. 

5. Conclusions 

The scintillation noise at Dome C for fast ex- 
posures is typically a factor 1.9 - 2.6 times lower 
than at Cerro Tololo and Cerro Pachon; leading to 
a corresponding reduction in this ultimate limit for 
high precision photometry. The "small aperture" 
scintillation index becomes important for adaptive 
optics when the distance between the actuators 
approaches the Fresnel zone size rp, and shadow 
patterns start becoming resolved (Masciadri et al. 
2004). Adaptive optics will also benefit from the 
long coherence time and large isoplanatic angle at 
Dome C, particularly in the infrared. 

For longer exposures, aj at Dome C is typi- 
cally 3.4 - 3.8 times less than at Cerro Tololo 
and Cerro Pachon. For a 60 second integration 
on a 4 m telescope, the median photometric er- 
ror is ^ 53 /zmag at Dome C. This parameter 
is important for exoplanet transit measurements 
because the change in fiux, caused by a transit- 
ing planet, is related to the planet Rp and star 
Rs radii by AF/F = (Rp/Rsf. For example, 
for a Jupiter size planet transiting a Sun size star 
AF/F = 0.01, for an Earth size planet this ratio 
is 0.0001. The lower scintillation noise at Dome C 
will allow for the transits of smaller planets to be 
detected than at the Chilean sites. 

The atmospheric contribution to the positional 
error in a differential astrometric measurement us- 
ing a long baseline interferometer at Dome C is al- 
ways less than the median values at Cerro Tololo 



Dome C Tololo 



Pachon 




200 400 600 800 1000 
Cj (arcsec rad"^m^/^s^/^) 





<r> 


6^ 


o 

T— ) 


>^ 




-(_) 




•'— ] 


in 




(^ 


X! 




crt 




X! 




n 




u 


o 


ft 


HJ 


(1) 




> 










lO 


(0 


C\2 


Hi 




a 




1=1 




CJ 


o 



L— — — ■ ■ ' ' 


^ — 


^- f~ - 


/ ■ '^ 






/ ■' ^ 




- 


/ ■ ^ 






/ ' ^ 






/ '^ 










- 


/ ..• 1 




1 



Cp (arcsec rad ' s ' ) 



Fig. 4. — Cumulative probabilities that the con- 
stants Ci and C2 are less than the value given, for 
Dome C, Cerro Tololo and Cerro Pachon. The 
astrometric error (Tatm (arcsecond) is equal to: 
Cii-i/20S-2/3 fo^ ^ ^ ^/Y and 9h <C B, and 

C2t-i/2^i/3 for t > eh/V and 9h > B. Here, t 
is the integration time, B is the baseline length, h 
and V are the turbulence-weighted effective atmo- 
spheric height and wind speed, and 9 is the stellar 
separation. 



Table 2 
isoplanatic angles and coherence times for adaptive optics 



Site 


9o(500 nm) to(500 nm) 


00 (2.2 Mm) 


To{2.2fim) 




arcscc 


ms 


arcsec 


ms 






Median values 






Dome C 


5.4 


7.2 


32 


42 


Cerro Tololo 


1.8 


2.0 


11 


12 


Cerro Pachon 


2.0 


2.6 

Average values 


12 


16 


Dome C 


5.9 


8.8 


35 


52 


Cerro Tololo 


1.9 


2.8 


11 


17 


Cerro Pachon 


2.1 


3.3 


13 


20 



Note. — tq and Oq at 500 nm are taken from the MASS data files, the 
values at 2.2 /im are scaled by A^/^. 
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Fig. 5. — The median error CTatm for three baseline 
lengths at Dome C, with an integration time of 1 h. 



and Cerro Pachon. This conclusion is obtained us- 
ing average wind profiles and remains provisional 
until a more complete analysis is done. 

Based on the expected low astrometric error at 
Dome C, a number of interferometric projects have 
already been proposed. These include the Antarc- 
tic Planet Interferometer (API; Swain et al. 2004), 
the Kiloparsec Explorer for Optical Planet Search 
(KEOPS; Vakih et al. 2004) and the Antarc- 
tic L-band Astrophysics Discovery Demonstrator 
for Interferometric Nulling (ALADDIN; Coude du 
Foresto et al. 2005). Many science programs would 
benefit from an Antarctic interferometer, includ- 
ing exoplanet detection and orbit determination 
and measurement of micro-lensing events (Lloyd 
et al. 2002). 
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thank James Lloyd and Mark Swain for partic- 
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of a "warm MASS" for making these measure- 
ments, and Victor Kornilov and Nicolai Shatsky 
for helpful advice on the data reduction software. 
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